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Measurement of D*"^ Cross Sections and the Charm Contribution to the 
Structure Function of the Proton in Deep Inelastic Scattering at HERA 

I. Redondo ^ *. On behalf of the ZEUS Collaboration 

''Dept. de Fisica Teorica, Universidad Autonoma de Madrid, Cantoblanco, Madrid 28049, SPAIN, 
e-mail: redondo@mail.desy.de 

We present measurements on D*^ production cross sections in deep inelastic e^p scattering with 1 < < 
600 GeV^ and 0.02 < y < 0.7 in two restricted kinematic regions in pt{D*^) and ri{D*^). The decay channels 
7?*"*" D^TT^ with subsequent decay K~-k^ or _R'~7r"'"7r+7r~ (+ c.c.) are used. The cross sections are 

extrapolated to the full kinematic region to determine the charm contribution to the proton structure function. 



1. INTRODUCTION 

The first HERA measurements of the charm 
contribution to the proton structure function, i^j*^ 
were reported by the HI and ZEUS Collabora- 
tions from an analysis of D*^ production in their 
1994 data sets [ |l|, ||. The results were con- 
sistent with Photon Gluon Fusion (PGP) being 
the dominant mechanism for D*^ production in 
e+p Deep Inelastic Scattering (DIS). If this is the 
case, this type of measurements are sensitive to 
the gluons in the proton. In addition, they can 
provide a test of the universality of the parton dis- 
tribution functions (pdf 's), namely, whether pdf 's 
extracted from the inclusive measurement of the 
proton structure function, F2, can be used as in- 
put for calculations of more exclusive processes as 
charm production. 

Here we present a study of D*^ production 
using the 1996 and 1997 data corresponding to 
an integrated luminosity of 37 pb^^. More 
than tenfold larger data sample, together with 
the modifications of the ZEUS detector made 
for the 1996 and 1997 operation allow an ex- 
tension of the kinematic range to both smaller 
and larger Q^. The D*^ is tagged via D*^ — s- 
(DO ^ A'-7r+)7r+(+c.c.) and D*+ 
A'~7r+7r+7r^)7r+(-|- c.c.) decay channels, refered 
to as Ar27r and Ar47r respectively. 
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2. D*^ CROSS SECTIONS 

The measured D*^ cross section using the 
Ar27r(Ar47r) final state, in the region 1.5(2.5) < 
piiD*) < 15 GeV, \r]{D*)\ < 1.5 is cr(e+p ^ 
e+D'^X) = 8M±OM{stat)1:oia{sys) nb(3.65± 
0.36{stat)t°oi°i{sys) nb). Figure |l| shows the 
differential D* cross sections in the restricted 
Q^, y, pt{D*) and r]{D*) region as functions 
of Zo5io(g2), logio{x), W, pt{D*), tj{D*) and 
x{D*)= 2p*{D*)/W, where p*iD*) is the mo- 
mentum in the 7*-proton CMS frame. The re- 
sults using each decay channel can be directly 
compared in the pt{D*) differential cross section. 
The agreement is satisfactory. 



3. COMPARISON WITH NLO QCD 

We compute NLO QCD calculations with a 
semi-inclusive Monte Carlo generator HVQDIS[ 
H for heavy quark production and subsequent 
fragmentation to D*^ via a Peterson fragmen- 
tation function [ ^ . This generator is based on 
NLO calculations [ ||] in the three flavor num- 
ber scheme (TENS), in which only light quarks 
(u, d, s) are included in the initial state pro- 
ton. Heavy quarks are produced exclusively by 
the convolution of the light flavours and the gluon 
with the massive matrix elements and coefficient 
functions calculated previously [ ||] . 

We use as input pdf ZEUS94 [ §. The QCD 
renormalization and factorization scales are set 
to ^yQ^+4m^ . rric is varied between 1.3 and 
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Figure 1. Differential cross sections for D* production from the K2tt final state (solid dots) and 
from the KAn channel (open triangles). The light shaded band corresponds to the standard Peterson 
fragmentation function. The dark shaded band show the NLO reweighted RAPGAP MC. The two sets 
of bands correspond to charm mass variations between 1.3 (upper curve) and 1.5 GeV (lower curve). In 
(a) and (b) the dark band is indistinguishable from the light one and is not shown. 



1.5 GeV. /(c D*) ^ 0.222 is taken from e+e~ 
measurements [ ^. The error on this quantity 
introduces a normalization uncertainty of ^ 9%. 
Finally, the Peterson fragmentation parameter is 
set to e = 0.035. The NLO QCD predictions (fig- 
ure |l|) are in reasonable agreement with the data, 
except in the Tj{D*) distribution, where the mea- 
surements show a shift into the positive 77 region 
(proton direction) with respect to the prediction. 
Also a softer charm fragmentation is favored by 
the data. 

3.1. Fragmentation Effects 

In Monte Carlo fragmentation models like JET- 
SET or HERWIG, a forward shift in the D* di- 
rection w.r.t. that of the original charm quark 
is produced during the fragmentation due to the 
interaction of the charm quark with the proton 
remnant via either strings or soft gluon radiatio 

To investigate how this affects the NLO QCD 
predictions we have reweighted a LO Monte Carlo 



for charm production, RAPGAP (it uses JET- 
SET for the fragmentation) in a way such that 
at the stage of the hard interaction it will repro- 
duce exactly the HVQDIS results for the pt{c) 
and 77(c) differential cross sections. The predic- 
tions from this NLO reweighted RAPGAP Monte 
Carlo are shown in (figure |l|). They provide a bet- 
ter descripion of the data, especially in the rj{D*) 
and x{D*) differential cross sections This result 
suggest that the small disagreement found with 
HVQDIS come from the fact that the Peterson 
function can not account for all the charm quark 
fragmentation effects present at HERA, in partic- 
ular the interaction with the remnant. 

4. Ff EXTRACTION 

The procedure for extracting starts by 
measuring the D*^ production cross section in 
the restricted pt{D*), 'ri{D*) region in bins of 
and y. The extrapolation to the full pt{D*), 
1]{D*) region is done in the following way: 
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with the ZEUS94 parametrization we determine 
that i^2^^ accounts for < 10% of F2 at low and 



Figure 2. The measured i^l"^ at constant as a 
function of The result of the NLO QCD ZEUS 
fit [ 01 is the solid curve. The dashed curves show 
the error from that fit. 



a\{e+p D*X) 



where Xi, is the center of gravity of the bin i 
and 'm' and 't' denote 'measured' and 'theoreti- 
cal' respectively. Ff* is taken from ZEUS94 [ § 
fit. For cr* we use the reweighted Monte Carlo. 
A number of assumptions are implicitly done in 
this procedure: 

• The TENS is valid, 

• F£'^ is negligible (< 1% of in our Q^, y 
region from calculations based on [ 

• the value of /(c — > D*) measured in e^e^ 
is valid also at HERA, 

• the cross section outside the restricted re- 
gion is well described by NLO QCD. 

Figure || shows the measured i^l'^ after combin- 
ing the results from both decay channels. Com- 
pared to our previous study we have extended 
the kinematic range to as low as 1.8 GeV^ 
and up to 130 GeV^ and the errors are reduced 
substantialy. J^|^ exhibits a steep rise with de- 
creasing a; at constant Q^. From a comparison 



for Q^> 11 GeV^ 



a; ~ 5 • 10--* and ~ 30% of F2 
at the lowest x measured. 

5. SUMMARY 

We have presented a charm analysis in DIS us- 
ing the combined 1996 and 1997 data sample. 
Charm was tagged with D* mesons decaying into 
two decays (K2t: and KAi:). In the experimen- 
tally accessible region, differential D*^ cross sec- 
tions are in reasonable agreement with NLO QCD 
calculations of charm production in the TENS us- 
ing a pdf extracted from an inclusive measure- 
ment of F2. This represents a successful test of 
the universality of the pdf 's. 

Small disagreements in the ri{D*) and x{D*) 
distributions show that the fragmentation a la 
Peterson can not account for all the charm quark 
fragmentation effects present at HERA. 

Using these calculations to extrapolate outside 
the measured pt(-D*), r]{D*) region, F^*^ was ex- 
tracted. Fl'^ is rising steeply with decreasing x at 
costant Q^. It amounts to ~ 25% of F2 at low x, 
Q2 > 10. 
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